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Abstract. The static dipole polarizabibility for the lsa electron state of the 
hydrogen molecular ion is calculated within Born-Oppcnhcimcr approximation. The 
variational expansion with randomly chosen exponents has been used for numerical 
studies. The results obtained for the dipole polarizability are accurate to the nine 
digits. 
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1. Introduction 

One of sample few-body systems, a H^" molecular ion have played vital roles in the 
development of molecular quantum mechanics. Moreover, this ion has been considered 
as a model system for the formulation of many different methods and approximations 
and is much easer to study theoretically than experimentally, and although its physical 
properties have been extensively calculated, few have been measured so far. 

At present, the system has been redeemed widely not only in theoretical field but 
in experimental sites. New type of variational expansions based on randomly chosen 
exponents [U [2] have been employed successfully by many researchers for the calculation 
of non-relativistic energy for the vibrational v = and rotational J = electronic 
state, even more the relativistic and radiative corrections of higher order of ma for 
lSa(v = 0, J = 0) state of the molecular ion [3j Hj. Within the Born-Oppenheimer 
approximation, the improved relativistic corrections of ma 4 and ma 6 orders for the 
ground state (v = 0, J = 0) of the H^" molecular ion have been obtained in [HI E]. 

There are, in addition, many theoretical precision nonadiabatic values of the 
electric-dipole polarizability of the ground state (v — 0, J — 0) of the molecular 
ion. Shertzer and Greene [Tj had evaluated the polarizabilities of the ground state 
(v — 0, J — 0) of the and using the finite element method. Bhatia and Drachman 
[8] calculated the dipole polarizability of the ground states (v = 0, J = 0) of the ions 
H.2 and using the second-order perturbation theory in Hylleraas coordinate. Taylor 
et al [9] obtained the electric-dipole polarizabilities of the ground state of H^" and D^" 
variationally using traditional molecular physics method in which any approximation 
based on the size of the electron mass relative to the nucleus mass could be avoided. 
Moss [I0] obtained rotationally averaged polarizabilities of the H^"(0,0) and D^(0,0) 
from scattering theory. Yan et al [TT] had derived more accurate theoretical values of 
the both nonrelativistic energy and dipole polarizability of ground state (v — 0, J = 0) of 
using variational wave function in Hylleraas coordinates. Hilico et al [12J had derived 
the most accurate results for the ground states (v = 0, J = 0) of hydrogen molecular ions 
H^" and using the perimetric coordinates. Korobov [13J had calculated the lowest 
order relativistic correction to the electric-dipole polarizability of the ground state of 
the molecular ion. 

Bishop and co-workers [HI [TH1 IS] had obtained the dipole polarizabilities of 
the ground state (v = 0, J = 0) of the ions in a clamped nucleus approximation. 
Various approaches for obtaining the polarizability in the adiabatic approximation were 
considered in [T7] . 

Static dipole polarizability of the lSa(v = 0, J = 0) electronic state of the the H^" 
molecular ion had been calculated theoretically many times in the Born-Oppeheimer 
approximation PS Q21 EQ1 ED E21 E3] • 

Furthermore, several laser spectroscopy experiments have been proposed [2^1 [251126] 
for high precision measurements of the rovibrational spectrum of hydrogen molecular 
ions and HD + . 
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The precision measurement of the scalar electric-dipole polarizability for the 
lSa(v = 0, J = 0) state of the molecular ion have been performed by Jacobson 
at al [271128]. 

In this work our goal is to obtain more accurate values for the static dipole 
polarizability for the lSa(v = 0, J = 0) state of than in [H] in the Born- 
Oppenheimer approximation. 

We want to show that the use of the variational expansion suggested in this 
work allows an analytical evaluation of the singular matrix elements required for the 
relativistic calculations and can provide us with very accurate data. 

In future studies it can be used to obtain the lowest order relativistic correction to 
the mean electronic dipole polarizability of the lSa(v = 0, J = 0) state of in the 
Born-Oppenheimer approximation. 

In Sect. 2, we outline the theoretical method for obtaining the the static dipole 
polarizability for the lSa(v = 0, J = 0) state of in the Born-Oppenheimer 
approximation. In Sect. 3 we consider the variational wave function and the 
nonrelativistic energy for this state of the H^~. Moreover, main results of numerical 
calculation will be presented in Tables in this section. In Sect. 4 we compare our 
results with the other values obtained in the Born-Oppenheimer approximation. Then 
a conclusion will be discussed. 

2. Theory 

The Schrodinger equation for the molecular ion may be solved at various different 
levels of approximation. In what follows, we will consider one of them; the Born- 
Oppenheimer approximation which means that an electron moves in the field of the 
clamped nuclei. 

In the clamped nucleus approximation, the electronic motion is perturbed by the 
electric field £, and the mean electronic dipole polarizability a(R) is derived as a function 
of the internuclear separation R. In this work the factor of (1 + e) = (2m + 2)/ (2m + 1) 
(details can be found in [TJElElQIllIIlQ^tlS], m is a mass of nucleus) in the perturbation 
will not be considered. Then a(R) is averaged over the rovibronic state. In the sum 
over rovibronic states approach, the adiabatic molecular wave function is written as 
a product of an electronic and nuclear function. In this paper, however, we do not 
consider the effects due to the rovibronic state, and think over only the mean electronic 
polarizability a(R) for ground state v — 0, J — of molecular ion. 

In the Born-Oppenheimer approximation the electronic wave function ty (r;R) 
satisfies the following Schrodinger equation 



where r*i and r^ are the distances from an electron to nuclei 1 and 2, respectively. 
Interaction with external electric field is given by 




(1) 



H' 



(2) 
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where d is the dipole moment and £ is electric field strength. 
Using the perturbation theory, it is easy to present that 

E = E + E"£ 2 (3) 

to the second order in £ , and E is the solution of Eq.([T]). 

If the electric field is parallel to the nuclear axis, the polarizability is expressed by 

a\\ = -2E" (4) 
and if it is perpendicular to the nuclear axis, the polarizability is expressed by 

a± = -2E" (5) 
The mean electronic polarizability is calculated by 

a = -(ay + 2aj_). (6) 
It can be shown [H, [29] that 

E" = BtDA (7) 
where B is the vector of coefficients in the expansion of the first-order wave function 

*' = £Wn (8) 

n 

(ip n is certain basic functions), A is the vector of coefficients in the expansion of the 
zero-order wave function 

n 

and D is the matrix elements with elements 

D nm = . (10) 

The B is obtained from the set of inhomogeneous equations: 

(H - E S)B = -DA (11) 
where the matrix elements of H and S are (ip' n Hoip' m ) and (ip' n il)' m ) , respectively 

3. Variational approximation and numerical results 

The variational wave function for the lsa(v = 0, J = 0) state of should be 
symmetrized and is constructed as follows 

oo 

$ = ai{e-^ ri -^ r2 + e- ftri ~ a ' r2 ). (12) 
i=i 

Parameters a, and are generated in a quasi-random manner 

a i =[U{i + l)^\{A 2 -A 1 )+A 1 (13) 

[x\ designates the fractional part of x, p a is a prime number, an interval [Ai,^] is a 
real variational interval, which has to be optimized. Parameters /?, are obtained in a 
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similar way. Details of the method and discussion of various aspects of its application 
can be found in [I] [2] . 

If the electric field is parallel to the nuclear axis, then 

cL = 



2R 1 

and the perturbed function ^' have a form 

oo 

% = X^( e ~ a ' ri_ftr2 + e~ Piri ~ aiT2 ). (14) 
i=i 

If the electric field is perpendicular to the nuclear axis, then 

d x = —r cos0 , 
and the perturbed function is taken the form 

oo 

tf' ± = rcos^^^e-" 1 " 1 -^ 2 + e - ftri - a * r ' 2 ) , (15) 

i=l 

where r is a distance from center along internuclear axis, 



r = — \j2r\r\ + 2r\R? + 2r\R 2 -r\-r\-R^. 

In Table 1 we show variational parameters employed in the calculation of energy values 
E and unperturbed wave function \l/ at internuclear distance R = 2. 0a.it.. Using this 
type of variational parameters leads to the very fast convergence. 



Table 1. Variational parameters and number of basic functions (ni) for different 
subsets of the variational wave function with N = 100. Intervals [A^Aa] and [B\, B2] 
correspond to part of a randomly chosen parameters on and [see Eq. (|13j) for dctials] , 
respectively, for the bond length R = 2.0 a. u. Prime numbers are p a = 2,pp = 3. 



N 


A 1 


A 2 


Bi 


B 2 


i=l 50 
i=2 50 


0.00 
0.20 


1.50 
2.00 


0.00 
1.00 


0.40 
6.00 



In Table 2 we present a comparison of the wave functions for the ground state 
(v — 0, J — 0) of the molecular ion at bond length R = 2.0 a.u.. The variationally 
obtained values of the \l/o wave function are more accurate than the previous ones. 

In Table 3, the convergence of both E (R = 2. 0a.it) and ^o(r = 2. 0a.it; R = 2. 0a.it) 
is presented. In this calculation the variational parameters presented in Table 1 are 
used. 

In Tables 4 the mean electronic polarizability for the lsa electronic state of 
molecular ion and comparison with [T3j are presented. The estimated accuracy of 
obtained values is nine significant digits. In Table 5 the convergence of an and a± 
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Table 2. The wave function for the lsa g electron state of the molecular ion at 
bond length R = 2.0 a. u. 



r 


Bates et al Gj 

7 E — 


— 3 


r 


Bates et al, 


30J 


*o(o-w.) 


0.00 


0.315 


0.314692 27 


1.60 






0.206 54247 


0.10 




0.315 970 96 


1.70 






0.180471 77 


0.20 




0.319 824 69 


1.80 






0.15761143 


0.30 




0.326 306 69 


1.90 






0.13758169 


0.40 




0.335 506 65 


2.00 


0.120 




0.120 04455 


0.50 




0.34755213 


2.10 






0.104699 83 


0.60 




0.362 610 66 


2.20 






9.128 154 54(-2) 


0.70 




0.380 892 46 


2.30 






7.955 43413(-2) 


0.80 




0.402 653 81 


2.40 






6.931 038 06(-2) 


0.90 




0.428 20125 


2.50 






6.036 634 71 (-2) 


1.00 


0.458 


0.457896 57 


2.60 






5.256 076 55(-2) 


1.10 




0.401762 60 


2.70 






4.575 156 24(-2) 


1.20 




0.352 20750 


2.80 






3.981 386 43(-2) 


1.30 




0.308 52410 


2.90 






3.463 80193(-2) 


1.40 




0.270 066 59 


3.00 


0.003 




3.012 782 38(-2) 


1.50 




0.236 248 99 


4.00 


0.007 




7.385 569 68(-3) 



Table 3. Convergence of both energy at bond length R = 2.0 a.u and wave function 
at distance from center along internuclear axis r = 2.0 a.u for the tsa g electron state 
of the H^" molecular ion. 

N E Q (a.u.) #o(a.u.) 

"60 -1.102 634 214 494 946 458 49 0.120 044 550 7422 
70 -1.102 634 214 494 946 460 60 0.120 044 550 805 9 
90 -1.102 634 214 494 946 46150 0.120 044 5508281 

100 -1.102 634 214 494 946 46150 0.120 044 550828 3 



with the increase in the expansion length of and ty' ± (Eqs. f|T4|) - f|T5|) ) is demonstrated 
for some values of internuclear distance R. 

Numerical evaluation of the matrix elements for operators in (j3J) and (JSD is 
expounded in the Appendix. 

In order to get accurate results we use three sets of basic function of the type (fT2j) 
(in a spirit of pQ) for small values of internuclear distance R, two sets for intermediate 
and large values of R, respectively. Total number of the basic function varies from 
N = 110 to N = 300. In our calculations arithmetics of sextuple precision (about 48 
decimal igits) implemented as a FORTRAN90 module has been used. In all tables the 
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Table 4. The dipole polarizability for the lsa(v = 0, J = 0) state of 
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factor x in the brackets means 10 x . Atomic units are used throughout. 
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Table 5. Convergence of dipole polarizabilities with the increase in the basic-set size 
of 



R(a.u.) 


N 


a\\ 


a± 


1.00 


100 


0.111959 391834 (1) 


0.776 229 394 339 




150 


0.111959 395159 (1) 


0.776 229 397334 




200 


0.111959 395159 (1) 


0.776 229 397 334 




250 


0.111959 395159 (1) 


0.776 229 397 335 




300 


0.111959 395159 (1) 


0.776 229 397 335 


1.50 


100 


0.244195 787807 (1) 


0.122 321061209(1) 




150 


0.244196 012 382 (1) 


0.122 321099 085(1) 




200 


0.244196 012 392 (1) 


0.122 321099 095(1) 




250 


0.244196 012 392 (1) 


0.122 321099 095(1) 




300 


0.244196 012 392 (1) 


0.122 321099 095(1) 


2.00 


100 


0.507 759 939 512 (1) 


0.175 763 635 406(1) 




150 


0.507 764 896 289 (1) 


0.175 764 860 893(1) 




200 


0.507 764 898 343 (1) 


0.175 764 861787(1) 




250 


0.507 764 898 353 (1) 


0.175 764 861793(1) 




300 


0.507 764 898 353 (1) 


0.175 764 861793(1) 



4. Conclusion 

The mean electronic polarizabibility for the lso~(v = 0, J = 0) electron state of the 
hydrogen molecular ion have been accurately calculated within Born-Oppenheimer 
approximation. The variational expansion with randomly chosen exponents has been 
used for numerical studies. This type of expansion allows us to use few number of 
basic functions. If nonrelativistic energy values of H^(0, 0) are required to 10~ 15 a.w. 
accuracy, then the wave functions must be accurate, at least, to this same level. The 
wave functions, however, are typically accurate to less than half as many significant 
figures as the energy. Seeing from Tables 3, the precision nonrelativistic energy value 
Eq at bond length R = 2.0 a.u can easily reach to 15-significant digits for basic function 
number N = 60 and the unperturbed wave function ty have 10-significant digits using 
the same number of terms. 

Previous calculations performed for lSo~(v = 0, J = 0) state of molecular ion 
over a wide range of internuclear separations R are by Dalgarno and Lewis [23], Calvert 
and Davison [21], Bates [2D] , McEachran and Smith |18j . and Bishop and Cheung 
[14] . Both Dalgarno and Lewis [23] and Bates [20j used the oscillator strength sum 
rule. McEachran and Smith [18] employed the variational procedure, using the accurate 
two-center James [31] orbital as the perturbed function. Bishop and Cheung p3] had 
calculated variationally the first accurate mean electronic polarizability over a wide 
range of internuclear separations. Rahman [19] and Adamov et al [22J had performed 
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Table 6. Comparison with earlier calculations at a bond length R — 2.0 a. u. 









pa 

E 1 


5.061 


1.758 


[23] 


5.283 


2.202 


m 


5.173 


1.847 


m 


5.084 


1.767 


L20J 


5.06 


1.76 


m 


5.199 


1.829 


m 


5.07765 


1.75765 


this work 


5.0776490 


1.7576486 



accurate first order variational calculations near the equilibrium distance R = 2.0 a.u.. 

In Table 6 we place the comparison of our results with the earlier ones, which 
demonstrate superiority of the newly obtained results. The results obtained for the 
dipole polarizability are accurate to the nine digits. That is two digits more accurate 
than in [14J. 

In general, the total molecular rovibronic wave function is taken a form as a 
product of separate electronic, vibrational and rotational wave functions, that is, 
^vJM = \v (J)) Yj(9, 0)-rovibrational wave function. Then the full static dipole 
polarizability can be written in the form 

a(v, J, M) = a e (v, J, M) + a v (v, J, M) + a r (v, J, M) (16) 

where the superscripts e, v, r imply the electronic, vibrational and rotational parts of 
the static dipole polarizability, respectively, and the electronic part-a e (t>, J,M) can be 
calculated by averaging over the mean electronic polarizability (details of Eq. (ITBTl can 
be found in [P7]). 

Seeing Eq. ( ITBl) . the rovibronic effects can be taken into account by averaging over 
the vibrational and rotational wave functions, of quantities calculated in the Born- 
Oppenheimer approximation. 

In the calculation of the mean electronic polarizability in this Letter the effects the 
electron mass respect to the nuclei mass has not been taken into account. This effects, 
however, must be taken into account when full three-body Hamiltonian (details cane be 
found [HO [32]) is employed. 

Using the unperturbed and perturbed ^' functions obtained in this paper we 
can calculate some expectation values of the lowest order relativistic correction (Eq. (9) 
in [13]) to the mean electronic polarizability a of the ISa state of molecular ion. In 
addition, the accurate data used in evaluation of the mean electronic polarizability a in 
this work can be employed in the calculation of relativistic correction to it. This work 
is in progress now. 
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Appendix: Analytical evaluation of the matrix elements 

The calculation of the matrix elements is reduced to evaluation of integrals of the type 
r lm (a,P)= [ r^r^e^-^dh. (A-l) 



Integers (I, m) are, in general, non-negative, but in case of singular matrix elements one 
of the indices can be negative. 

The function r 00 can be easily obtained 

47r e~P R — e~ aR 

r 00 (a,(3,R) = — a2 _ p2 , (A-2) 

where R is the distance between nuclei, then T [m (a, P; R) for non-negative (l,m) may 
be generated from ( 1A-2I) by means of relation 

d V ( d xm 



r^;*)=(--l (--1 T 00 (a,f3,R). (A-3) 
Integral r_i («,/3; R) is expressed by 

T_ 10 (a,(3;R) = ^{e" fl E 1 (i?(a + (3)) + e^ R \n R(a + P) 

RP L (A-4) 

-e" fl [Ei(i2(a - P)) + \nR(a - /?)]}. 

Worthy to note that a function in square brackets is analytic when argument is 
zero. Integrals r_ l m are generated from r_i j0 similar to (1A-3|) : 

r_ 1)m (a, p ] R)=(-—) r_ 1)0 (a, /5, i?). (A-5) 



9/? 

Function Ei(^) encountered in (|A-4[) is the exponential integral function |33j : 



/»oo 

Ei(z) = T(0,z) = / r x e- l dt. 

J z 
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